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capacity  distributions  with  respect  to  two  physically  important  distances,  the  depth  and 
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Scientific  Progress 

Computer  Science,  Marshall  Lochbaum  and  Bihan  Zhang,  undertook  to  construct  from 
plant  root  images  a  3D  piecewise  smooth  model  for  the  plant  root  system  with  the 
additional  data  of  varying  cross-sectional  area  of  the  roots.  The  resulting  geometric 
and  physical  properties  of  the  plant  root  system  would  then  be  compared,  using 
programs  which  the  investigator  had  previously  written,  with  several  theoretical  models 
involving  the  intake  capacity  distributions  with  respect  to  depth  and  the  radial  distance 
from  the  vertical  center  line.  In  particular,  certain  properties  the  investigator  had 
observed  in  several  cases  he  had  analyzed  by  hand,  would  be  compared  with  results 
for  a  time  series  for  a  single  plant  root  system,  and  further  compared  with  results  for 
multiple  plants  of  the  same  species  and  across  a  number  of  species. 

The  construction  of  the  3D  model  is  derived  from  a  collection  of  20  images  at 
successive  18  degree  angles  about  a  central  vertical  axis.  The  approach  to  construct 
the  model  was  a  two  step  process  where  first  a  3D  volumetric  model  would  be 
constructed,  and  then  a  piecewise  smooth  model  with  cross-sectional  data  would  be 
constructed  from  the  volumetric  model. 

The  investigator  consulted  with  Gary  Bishop  and  Leonard  McMillan  who  both 
contributed  to  the  development  of  the  “visual  hull  method”  for  reconstructing  a  3D  object 
from  multiple  projections.  They  pointed  him  to  a  library  of  software  by  Matusik  for 
implementing  the  method.  After  performing  several  alignment  procedures  to  the 
images,  it  was  discovered  that  experimental  method  for  capturing  the  images  did  not 
include  accurate  “camera  calibration”  which  is  a  basic  ingredient  for  the  most  effective 
use  of  the  software.  This  forced  a  basic  reevaluation  of  how  to  apply  the  method  and 
compensate  for  this  gap.  Also,  the  visual  hull  method  is  especially  effective  for  objects 
that  are  largely  convex;  however,  the  many  roots  in  the  plant  root  image  produced  a 
considerable  amount  of  noise,  which  made  the  resulting  volumetric  model  of  limited  use. 
Much  effort  went  into  trying  methods  to  overcome  these  problems.  These 
included  attempts  to  identify  from  among  the  images  a  subset  which  minimized  noise. 

A  second  method  was  applied  which  employed  a  Markov  Chain  Monte  Carlo  method  for 
sampling  the  3D  viewspace  for  points  in  multiple  images  to  construct  a  probabilistic 
distribution  for  the  3D  volumetric  model.  This  was  to  be  combined  with  a  section 
method  for  identifying  at  identified  root  points  an  elliptical  cross-sectional  model,  which 
together  with  root  direction  would  give  an  orthogonal  cross-section. 

Various  parts  of  this  process  were  implemented  in  code  written  in  either  Matlab, 
or  the  language  “J”.  The  work  still  had  not  reached  the  point  where  it  was  possible  to 
identify  all  of  the  roots  in  the  image  without  introducing  by  the  user  additional  seed 
points  for  the  method.  Also,  only  preliminary  steps  were  taken  to  move  from  the 
probabilistic  model  of  point  densities  to  obtain  an  accurate  volumetric  model.  There 
then  remained  the  issue  of  constructing  the  piecewise  linear  model  central  to  the 
volumetric  one  and  using  the  elliptical  cross-section  to  complete  the  model. 

Time  became  an  important  limiting  factor  in  the  project  because  of  the  delay  in 
obtaining  the  computer  on  which  was  loaded  the  images,  compilers  and  software  which 
we  used.  The  grant  began  at  the  end  of  October,  with  the  computer  being  ordered  then; 
however,  the  Purchasing  Dept,  at  UNO  did  not  deliver  the  computer  until  early 
December,  at  which  time  the  undergraduates  were  starting  their  final  exams,  followed 
by  the  Christmas  vacation.  This  meant  that  except  for  some  preliminary  discussions, 
the  work  did  not  really  begin  until  January,  and  continued  until  roughly  mid  May,  when 
the  academic  year  ended  and  the  students  left  for  summer  jobs. 

The  investigator  had  hoped  that  the  project  would  have  been  further  enough 
along  by  June  so  that  he  could  continue  using  what  had  been  developed  to  begin 
analyzing  the  variety  of  root  images  and  testing  the  results  against  the  theoretical  ideas. 
Because  this  was  not  possible,  the  investigator  asked  that  almost  all  of  the  funds  from 
the  grant  (  $44,793  of  the  original  grant  award  $50,000)  be  de-obligated  It  was  intended 
to  reapply  for  a  STIR  grant  for  the  next  academic  year  requesting  the  remaining  funds  to 
continue  with  the  project.  Because  the  investigator  received  a  University  leave  to  be 
away  the  next  academic  year,  the  reapplication  was  postponed,  but  has  not  been 
abandoned. 
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Capacity  Distributions  for  Plant  Root  Structures” 

In  this  “STIR”  grant,  this  investigator  with  two  undergraduate  students  in 
Computer  Science,  Marshall  Lochbaum  and  Bihan  Zhang,  undertook  to  construct  from 
plant  root  images  a  3D  piecewise  smooth  model  for  the  plant  root  system  with  the 
additional  data  of  varying  cross-sectional  area  of  the  roots.  The  resulting  geometric 
and  physical  properties  of  the  plant  root  system  would  then  be  compared,  using 
programs  which  the  investigator  had  previously  written,  with  several  theoretical  models 
involving  the  intake  capacity  distributions  with  respect  to  depth  and  the  radial  distance 
from  the  vertical  center  line.  In  particular,  certain  properties  the  investigator  had 
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The  construction  of  the  3D  model  is  derived  from  a  collection  of  20  images  at 
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the  model  was  a  two  step  process  where  first  a  3D  volumetric  model  would  be 
constructed,  and  then  a  piecewise  smooth  model  with  cross-sectional  data  would  be 
constructed  from  the  volumetric  model. 

The  investigator  consulted  with  Gary  Bishop  and  Leonard  McMillan  who  both 
contributed  to  the  development  of  the  “visual  hull  method”  for  reconstructing  a  3D  object 
from  multiple  projections.  They  pointed  him  to  a  library  of  software  by  Matusik  for 
implementing  the  method.  After  performing  several  alignment  procedures  to  the 
images,  it  was  discovered  that  experimental  method  for  capturing  the  images  did  not 
include  accurate  “camera  calibration”  which  is  a  basic  ingredient  for  the  most  effective 
use  of  the  software.  This  forced  a  basic  reevaluation  of  how  to  apply  the  method  and 
compensate  for  this  gap.  Also,  the  visual  hull  method  is  especially  effective  for  objects 
that  are  largely  convex;  however,  the  many  roots  in  the  plant  root  image  produced  a 
considerable  amount  of  noise,  which  made  the  resulting  volumetric  model  of  limited  use. 

Much  effort  went  into  trying  methods  to  overcome  these  problems.  These 
included  attempts  to  identify  from  among  the  images  a  subset  which  minimized  noise. 

A  second  method  was  applied  which  employed  a  Markov  Chain  Monte  Carlo  method  for 
sampling  the  3D  viewspace  for  points  in  multiple  images  to  construct  a  probabilistic 
distribution  for  the  3D  volumetric  model.  This  was  to  be  combined  with  a  section 
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Various  parts  of  this  process  were  implemented  in  code  written  in  either  Matlab, 
or  the  language  “J”.  The  work  still  had  not  reached  the  point  where  it  was  possible  to 
identify  all  of  the  roots  in  the  image  without  introducing  by  the  user  additional  seed 
points  for  the  method.  Also,  only  preliminary  steps  were  taken  to  move  from  the 
probabilistic  model  of  point  densities  to  obtain  an  accurate  volumetric  model.  There 
then  remained  the  issue  of  constructing  the  piecewise  linear  model  central  to  the 
volumetric  one  and  using  the  elliptical  cross-section  to  complete  the  model. 


Time  became  an  important  limiting  factor  in  the  project  because  of  the  delay  in 
obtaining  the  computer  on  which  was  loaded  the  images,  compilers  and  software  which 
we  used.  The  grant  began  at  the  end  of  October,  with  the  computer  being  ordered  then; 
however,  the  Purchasing  Dept,  at  UNO  did  not  deliver  the  computer  until  early 
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